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Abstract. Using the Belle detector we study the characteristics of beam collisions 
at the KEKB 3.5GeVe+ on 8GeVe~ asymmetric energy collider. We investigate the 
collision timing tip and its z-coordinate along the beam axis zjp as a function of the 
position of the colliding bunch in a beam train. The various tip and zip behaviors 
observed by Belle are attributed to beam loading effects in the radio frequency cavities 
that accelerate the beams with a beam abort gap. We report these results in detail 
and discuss the prospects for the SuperKEKB collider. 
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1. Introduction 

The B-factories, Belle [I] at KEKB [2] and BaBar [3] at PEPII [1], have carried out 
precision measurements of various decay modes of B mesons to verify the Kobayashi- 
Maskawa mechanism [5] that explains CP-violation in the Standard Model. Belle is a 
sophisticated Air detector, which provides precision measurements of the S-meson decay 
vertex, excellent momentum resolution and particle identification capability. KEKB is 
an asymmetric collider with beams of S.SGeVe"'' in the low energy ring (LER) and 
8GeVe~ in the high energy ring (HER). The BaBar detector at the PEPII collider is 
comparable to the Belle detector at the KEKB collider. The total integrated luminosity 
of the two B-factories is 1.5 ab~^ after ten years of data taking. 

KEK plans to upgrade the detector and collider to study CP-violation in 5-decays 
with extremely high precision to search for new physics beyond the Standard Model, 
with an integrated luminosity of at least 50ab~^ [6j. The upgrade plan is to develop 
an asymmetric e^e~ collider with a super high luminosity of 8 x 10^^cm~^s~^, and to 
develop a detector with improved performance under higher beam backgrounds (due to 
beam currents that are a few times larger). There are significant challenges to developing 
such a SuperKEKB factory [7j. 

We report a study of the behaviors of the beam collision timing tjp and its position 
Zip along the beam axis as measured by Belle. These behaviors are finally attributed to 
the beam phases of the colliding bunches in KEKB. To date, there has been no report 
on these behaviors, as it requires a collaboration between experts from the detector and 
the collider. Moreover, the accelerator parameters relevant in this study were measured 
only during certain periods of machine study, while the Belle measurement is available 
throughout physics data taking. In this study, we use the time-of- flight detector (TOF), 
which has a time resolution of 110 ps, for the study of the collision timing tip and the 
silicon vertex detector (SVD) and the central drift chamber (CDC) for the study of 
the collision 2;-coordinate zip, which is measured with a resolution of about 100 /xm [8]. 
These Belle measurements of tip and zip can be understood in terms of KEKB design 
parameters and operating conditions. 

Our goal is to understand the sources of systematic uncertainties in the timing 
measurements with Belle at KEKB, and to estimate prospects for the Super KEKB- 
factory, in particular for particle identification detectors such as the time-of-propagation 
(TOP) [HI |9] and Cherenkov timing detectors [HI HO], which require precision timing 
measurements with systematic uncertainty better than 20 ps. 

1.1. Belle detector 

A comprehensive description of the Belle detector is given elsewhere [1]. Figure [1] 
shows the Belle detector configuration. Charged tracks are reconstructed from hit 
information in a central drift chamber (CDC) located in a uniform 1.5 T solenoidal 
magnetic field. The z axis of the detector and the solenoid are along the positron 
beam, with the positrons moving in the negative z direction. The CDC measures the 
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longitudinal and transverse momentum components (along the z axis and in the r0 
plane, respectively). Charged particle trajectory coordinates near the collision point are 
provided by the SVD. Photon detection and energy measurements are performed with 
a CsI(Tl) electromagnetic calorimeter (ECL). The ECL is also used for identification 
of electrons. Identification of pions, kaons and protons is performed using information 
from the time-of-flight counters (TOP) and silica-aerogel Cherenkov counters (ACC). 
The ACC provides good separation between kaons and pions or muons at momenta above 
1.2 GeV/c. The TOP system consists of a barrel of 128 plastic scintillation counters and 
can distinguish K/ir for tracks with momenta below 1.2 GeV/c. Lower energy kaons are 
also identified using specific ionization (dE/dx) measurements in the CDC. The magnet 
return yoke is instrumented to form the and muon detector (KLM), which detects 
muon tracks and provides trigger signals. 

In order to study the characteristics of tip and zip, we use a large sample of /i- 
pair events (e"^e~ — )■ which were recorded using the standard Belle trigger and 

selected through the standard Belle data processing. The accumulation rate of /i-pair 
events is about 9 Hz at a nominal luminosity of 1.5 x 10^^cm~^s~^. The purity of the 
selected //-pair events is better than 98%. 




Figure 1. The Belle detector. 

1.2. KERB 

The KEKB collider achieved a peak luminosity of 2.1 x 10'^'^cm~^s~^ in June 2009, which 
is more than twice the design value of 1 x 10^'^cm~^s~^. An integrated luminosity of 
1018 fb~^ was accumulated with the Belle detector by December 2009. 
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Figure |2] shows the KEKB colhder complex [2], which has a circumference of 
approximately 3 km. KEKB stores beam currents up to 2.0 A in the LER and 1.35 A 
in the HER with 508.9 MHz radio frequency (RF) acceleration systems, which are 
located at the FUJI, OHO and NIKKO straight sections. The RF systems consist 
of two types of higher order modes damped cavities: the normal-conducting Accelerator 
Resonantly-coupled with Energy-Storage (ARES) cavities, and Superconducting (SC) 
cavities. Twenty ARES cavities are used in the LER to handle higher beam currents, 
while twelve ARES and eight SC cavities are used in the HER to achieve a higher 
voltage [11]. The RF reference signal, generated at the KEKB control room near the 
FUJI hall, is distributed through two coaxial cable lines in the clock- and anti-clock-wise 
directions around the KEKB collider. These cables have good phase versus temperature 
characteristics, less than 10 ppm/°C. The RF reference signals are phase-locked at each 
station and distributed to the RF cavities. Thus, the accuracy of the 508.9 MHz RF 
reference clock is found to be typically 6 ps in one circulation around the KEKB collider. 
The Belle detector uses the same RF reference signal distributed at TSUKUBA hall. 

The LER and HER beams collide bunch-by-bunch at the IP with a finite angle of 
22 mrad in the horizontal plane. Each beam train consists of 5120 buckets synchronous 
to the 508.9 MHz RF clock. In usual operation, the buckets are filled with three- or 
four-bucket spacing and with a period of 49 RF buckets. A typical beam train contains 
1584 bunches (buckets filled by beam). The beam train is followed by a beam abort gap 
(BAG) of 270 empty buckets, which is needed to safely kick out the beam from the ring. 
Figure [3] shows a typical collision bunch pattern as a function of bucket number, which 
is obtained from analysis of /i-pair events accumulated with Belle. Here the number of 
/i-pair events is shown as a function of the bucket number. We clearly see the bunch 
train followed by the BAG as described above. 

In December 2006, two superconducting crab cavities [12] were installed in the 
LER and the HER, one for each, in order to operate KEKB in a crab crossing scheme 
aiming to increase luminosity. Figure |4] illustrates (a) finite angle collisions and (b) crab 
crossing at the IP. In the crab crossing scheme the bunches are kicked in the horizontal 
plane by transverse RF in the crab cavities (located in the NIKKO straight sections) so 
that they rotate and collide head-on at the IP (TSUKUBA) even in finite-angle crossing 
beam orbits. 

1.3. Belle data sets used in this analysis 

We study three data sets that are summarized in Table [H by the run period, LER and 
HER beam currents, recorded peak luminosity, CM energy and status of crab cavities. 

(1) Exp55 is the data set before crab cavity installation (no-crab crossing), in which 
KEKB achieved the highest peak luminosity with the highest LER and HER currents. 

(2) Exp61 is the data after the crab cavity installation (crab crossing), during which 
KEKB was operated with head-on collisions. (3) Exp65 is the Spring 2008 run period, 
in which we recorded data at the center-of-mass (CM) energies of the T(4S) and T(1S) 



KERB Beam Collision Stability observed with the Belle Detector 



5 




Figure 2. The KEKB collider eomplex. 




Figure 3. The number of /i-pair events as a function of the RF bucket number, which 
are reconstructed by the Belle detector for Exp61. See Appendix [Appendix A.l| for 
details. 



resonances. To change the CM energy, both the LER and HER beam energies are 
changed proportionally. In Table [H T(1S), T(4S) and T(5S) indicate the CM energies 
of 9.460 GeV, 10.580 GeV, 10.860 GeV, respectively. Here T(6S) stands for the energy 
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HER LER HER LER 

(a) Finite angle collision (b) Crab collision 



Figure 4. Collision schemes at KEKB: (a) finite angle collisions and (b) crab crossing, 
region just above the T(5S). We use more than 300k /i-tracks for each data sample. 

Table 1. Summary of the three run periods (Exp55, Exp61 and Exp65) used in this 
analysis. The peak luminosity is in units of 10'^^ cm^^s^^. 





Run period 


LER 


HER 


Peak lum. 


CM energy 


Crab 


Exp55 


Sept.-Dec. 2006 


1.65A 


1.37A 


1.7 


T(4S) 




Exp61 


Oct.-Dec. 2007 


1.62 A 


0.81 A 


1.4 


T(4S,5S,6S) 


ON 


Exp65 


Apr.- June 2008 


1.61 A 


0.95A 


1.6 


T(4S,1S) 


ON 



2. Collision timing and collision position 

2.1. Variation of Collision timing tjp 

We first discuss the long term behavior of the collision timing tip. We define tip=(tobs — 
texp) using /x-pair events. Here, tobs is the measured time-of-fiight for one fj, candidate and 
texp is the prediction from its track trajectory in the CDC (see Appendices Appendix A.l 



and [Appendix A.2| ). We usually obtain a mean tip=tobs — ^cxp = for each run using 
the standard TOP calibration after run-by-run correction. Por online data processing, 
we use a default database obtained by TOP calibration using /i-pair events from the 
previous run period. Once the database is installed at the beginning of the new run 
period, it is unchanged during the whole run period, typically for several months. If 
both the KEKB collider and the Belle detector (TOP and CDC) are stable, the mean 
tip value should be unchanged. We conclude that the variation of the online tip value 
must be due to changes in beam conditions or changes in the RP reference timing from 
KEKB to the TOP readout system. 

Pigure [5] shows the mean tip as a function of run number for Exp65, taken at 
the energies of the (a) T(4S) and (b) T(1S) from April to June 2008. We see a slow 
drift within a ±10 ps range and step-like shifts at runs 610 and 810. The tip changes 
are due to various effects in the KEKB collider and the Belle detectors. While the 
individual sources are not understood completely, the tip variation can be corrected 
even with significant changes of beam energy and frequent optics corrections during 
KEKB operation. In practice, we apply a run-by-run correction for tip to reduce the 
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bias for physics analysis. We have a time resolution of 110 ps for the TOF system 
(including a readout system error of 40 ps) after ten years of operation, starting with a 
value of 96 ps in 1999 (see Appendix Appendix A.2). 




2GG 400 600 800 1050 1100 1150 1200 

Run Number (Exp65) Run Number (Exp65) 



Figure 5. tip versus run number for Exp65. (a) T(4S) runs and (b) T(1S) runs. The 
run number is reset to 1001 for T(1S) runs after run 810. 

In the analysis that follows, we use the tip value averaged over a run period and 
corrected for the run-by-run variation. Here, a run typically extends for a few hours to 
a day and contains about 300 k /x-pair events. We normally assume that the collision 
timing is independent of the bunch position in a beam train. However, we observe a 
fine structure in the collision timing as discussed in the following sections. 

2.2. tip and zjp vs bucket number by Belle 

The position of each beam bunch may be displaced from that synchronized with the 
RF frequency, because of dynamic beam loading effects in the RF cavities and the 
longitudinal impedance of other accelerating components. Therefore, it is interesting to 
study tip and zip as a function of the bucket number A^bucket, the collision bunch position 
in a bunch train. 

Figure [6] shows the Belle measurements of tip and 2;ip as a function of Abuckct; (a) 
and (c) show the results for Exp55 while (b) and (d) are the corresponding results for 
Exp61. The open circles indicate the Belle measurements. We see linear tip variation 
both for Exp55 and Exp61. On the other hand, we observe parabolic zip behavior 
for Exp61, while zip is almost flat for Exp55. These Belle observations of tip and 2;ip 
behaviors are discussed quantitatively in the next section. 
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Figure 6. tip and zjp vs bucket number; (a) and (c) for Exp55, and (b) and (d) for 
Exp61, respectively. The open circles are Belle measurements, while the filled circles 
are predictions from the measured LER and HER beam phases (see text). 



2.3. Beam phase and transient beam loading at KERB 

Figure [7] shows the beam phases as a function of the bucket number (Nbucket) measured 
by KEKB in the same periods of (a) Exp55 and (b) Exp61. The beam phase is 
defined as the phase with respect to the corresponding RF bucket. Here, the phases are 
measured at every 49-th bucket by the gated beam position monitor (GBPM) [13] (see 



Appendix Appendix B.l ), and the phase of the first bunch after the beam abort gap 
(bucket number 0) is set to zero by definition. We observe a beam phase advance with 
increasing bucket number in the LER and the HER with a maximum phase change of 
approximately 3.5 degrees. Figure [7] (a) shows a similar pattern of beam phases in the 
LER and the HER for Exp55, while Fig. [7] (b) shows a systematic difference between 
them for Exp61. 
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Figure 7. Beam phase as a function of bucket number for (a) Exp55 and (b) Exp61. 



The observed beam phase modulation is known as "transient beam loading" in the 
RF cavities due to the existence of a beam abort gap (BAG) after a bunch train. The 
beam phase change between the head and tail bunches is approximately proportional to 
the total beam current and the length of the BAG. The modulation pattern is sensitive to 
the type of RF cavities (ARES or SC) and the total beam current, as well as the pattern 



of bunch currents in the bunch train (see Section [3] and Appendix Appendix B.2) 



2.4- Effect of beam phases to tip and zip 

In this section, we give a quantitative discussion of the tip and zjp behaviors observed by 
Belle from beam phases of the collision bunches in the beam train. Assuming head-on 
collisions of the LER and HER bunches in the horizontal plane at the interaction point 
(IP) for a finite-angle collision with crab crossing, we derive the following formulae, 

' ^HER + 0LEr) 1 



tip — Trf X 



Zip — Lrf X 



^HER 



^LERj 



360 
1 



ns 



m 



(1) 



(2) 



2 360 

where Trf (approximately 2 ns) is one RF period and Lrf (approximately 0.6 m) is one 
bucket interval given by cTrf where c is the speed of light. The beam phases 0her 
and 0LER are in degrees. Here we assume that the change of zjp is much smaller than 
the vertical beta function /3* of 6 mm at the IP. If this is not the case, the change of 
Zip will be smaller than the value given by Eqn. (|2]) due to the hourglass effect. The 
collision timing of the tail bunch is advanced by about 20 ps relative to the synchronous 
RF bucket, so that tip is delayed by the same amount. It should be noted that the zip 
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behavior is very sensitive to the LER and HER beam phase difference, while the tip 
variation is almost proportional to the beam phase modulation. 

In finite-angle collisions without crab cavities, the bunches are more separated in the 
horizontal plane before and after the collision than in the case of the head-on collisions 
with crab cavities, as seen in Fig. HI We estimate the zjp variation by a naive calculation 
from the geometrical overlap of the colliding LER and HER bunches in three dimensions 
as a function of time, and taking into account finite-angle collisions with and without 
crab cavities and the beam focusing optics at the IP. We find that the zip variation is 
slightly smaller, about 70% of that predicted by Eqn. (|2]). 

Figure [6] shows the Belle measurements (open circles) of tip and zjp compared with 
the predictions (filled circles) from Eqns. ([1]) and ([2]). They are in good agreement except 
for a small systematic difference in Fig. [6](d). In Exp55, the LER and HER beam phase 
variations are similar, so that the zip behavior is flat and has no distinctive structure. 
In the case of Exp61, we see a parabolic zip behavior with a range of approximately 
0.6 mm, due to the significant difference between the LER and HER beam phases. This 
is attributed to the difference in the LER beam phase behaviors between Exp55 and 
Exp61 as shown in Fig. [71 The cause of the difference is discussed in Section [3l 

Thus, the observed tip and zip behaviors are explained well by the beam phase 
change as a function of the bunch position in the colliding bunch trains (beam phase 
modulation), which is caused by the transient beam loading effect in the RF system due 
to the beam abort gap in the beam train. 

2. 5. Size of Collision point zpp 

Figures [8] (a) and (b) show the zip distributions for Exp55 and Exp61, respectively. 
The widths of these distributions, fitted to a Gaussian function with a = 3.71mm and 
3.96 mm, respectively, are affected by the LER and HER bunch sizes along the beam 
axis and the collision scheme (no-crab or crab crossing). 

The natural bunch lengths for (LER, HER) are estimated to be (4.86 mm, 4.84 mm) 
for Exp55 and (4.86 mm, 5.22 mm) for Exp61 from the known KEKB operation 
parameters (Table [2]). The last two columns give the expected zip widths and the 
corresponding Belle measurements. Here, the zip width is estimated from the natural 
bunch lengths by crip = 1/ \/l / o"ler + V "^her- The widths measured by Belle are 
about 10 % larger than these estimates, due to the bunch lengthening effect caused by 
the inductive part of the longitudinal impedance of the ring components. The zip widths 
in the crab crossing runs (Exp61) seem to be larger than those in the no-crab crossing 
runs (Exp55). 

3. Discussion on the LER beam phase behaviors 

Here, we discuss the cause of the significant difference in the LER beam phase 
behaviors in Exp55 and Exp61 due to beam loading in the RF cavities (see 
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Figure 8. zip distributions for (a) Exp55 and (b) Exp61. The Gaussian widths are 
3.71 mm and 3.96 mm, respectively. 



Table 2. Size of coUision position zjp and LER and HER bunch sizes 





natural bunch length 


Zip width 


zjp width 


crab 
status 


LER 


HER 


KEKB est. 


Belle meas. 


Exp55 


4.86 mm 


4.84 mm 


3.43 mm 


3.71 mm 


OFF 


Exp61 


4.86 mm 


5.22 mm 


3.56 mm 


3.96 mm 


ON 



Appendix Appendix B.2 ). We first consider the possibility that the striking difference is 
related to the crab cavity installation. However, as crab cavities of the same type were 
installed in both the LER and HER, we cannot explain why we observed a significant 
difference only in the LER beam as shown in Fig. [3 

In general, the beam phase change with bunch position is attributed to transient 
beam loading in the RF cavities, which is sensitive to the total beam current and the 
bunch current pattern after the BAG. The transient beam loading effect is characterized 
by the ratio Tq/Tj, where Tq is the beam revolution period and Tf is the filling time of 
the RF cavity. In the limit Tq S> Ty, the beam phase shift (relative to the RF bucket) is 
expected to be largest in the beam head just after the BAG. This results in non-linear 
behavior of the beam phase. On the other hand, in the limit Tq <^ Tf, the beam phase 
is expected to be linear with bunch position as the effect is averaged over many beam 
revolutions. ARES cavities have a Tq/Tj value of 0.54, which is about twice as large as 
SC cavities. Therefore, the LER, which consists of ARES cavities only, tends to have 
stronger non-linear behavior of the beam phase. 

The above discussion is based on the assumption that the bunch current (current 
intensity) is uniform along the beam train. Figure [9] shows the LER and HER bunch 
current variation versus bunch position for (a) Exp55 and (b) Exp61. We see that the 
HER bunch currents are constant in both; the HER beam phases vary almost linearly 
with bunch position as shown in Fig. [71 Using Eqn. f lB.2p . we calculate the ratio of the 
head to tail phase shifts to be 0.74 for Exp61/Exp55 from an HER current of 0.80 A 
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Figure 9. Bunch current at each 49-th bucket for (a) Exp55 and (b) Exp61. 



(1.26 A) and an RF voltage of 13 MV (15 MV). We also measure an HER beam phase 
shift of 2.8 (Exp61) and 3.7 (Exp55) degrees, and find the ratio of 0.76, which is in good 
agreement with the above estimate. 

On the other hand, the LER bunch current in Exp55 shows a slow increase of about 
20% in the head part (up to bucket number 1000), while it is almost uniform in Exp61. 
We conclude that the non-linear behavior of the beam phase (Exp61) can be attributed 
to the transient response of the ARES cavities to the uniform bunch currents along the 
bunch train. In the case of Exp55, the pattern of non-uniform bunch currents weakens 
the steep increase of the beam phase in the head part and increases the phase change 
in the tail part, and thus gives almost linear phase behavior. 

The SuperKEKB [7] collider was initially designed to store currents of 9.4 A in the 
LER and 4.1 A in the HER (High Current scheme), which are four times larger than 
those of KEKB. To reduce the beam-loading effect due to the beam abort gap, which is 
proportional to the beam currents, we have designed a gap length of 2% (5% at KEKB) 
that will require development of new beam abort kickers. Recently an alternative design 
called the Nano-Beam scheme has been proposed and is the baseline option. This scheme 
obtains high luminosity via small beam sizes, which are achieved by reducing emittances 
and /3-functions. The required beam currents of 3.3 A (LER) and 1.9 A (HER) are 
smaller than in the initial SuperKEKB proposal and hence the beam loading effects are 
somewhat reduced. 

4. Summary 

Using the Belle detector, we study the characteristics of colliding beams at the KEKB 
asymmetric energy collider. The behaviors of the collision timing (tip) and its z- 
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coordinate along the beam axis (-^ip) observed by Belle are well explained by the beam 
phase variations along the LER and HER beam trains (beam phase modulation). Here, 
we calibrate for the long term variations of tip (see subsection 12.11) whose causes are 
not well understood. The studied beam phase variations are caused by transient beam 
loading in the RF systems due to the beam abort gap in the beam train. The zip 
behavior is very sensitive to the beam phase difference between LER and HER beam 
trains. The beam phase variation is sensitive to the bunch current pattern along the 
beam train as well. In collisions of LER 1.6 A and HER 1.2 A beams with a 5% beam 
abort gap at KEKB, we find an Zip variation of about 0.6 mm, and an tip variation of 
about 20 ps. To reduce the beam loading effect, we have designed a gap length of 2% 
to allow the SuperKEKB collider to store currents of 9.4 A (LER) and 4.1 A (HER) for 
the High Beam current scheme. The beam loading effect is much reduced with the same 
gap design for the Nano Beam scheme. 

We are able to reduce the systematic errors in collision timing and position due to 
the beam loading in the RF system to insignificant levels by correcting the data for the 
fine structure as a function of the bucket number, as discussed in this report. A crucial 
issue for the SuperKEKB factory detector, which must realize a systematic error of less 
than 20 ps, will be to reduce time jitter in the readout system synchronized with the 
RF clock and in the calculation of the time-of-flight value from tracking information. 
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Appendix A. Event timing determination with the TOF detector 

Appendix A.l. Bucket number determination 

A detailed description of the TOF readout system and analysis procedure is found 
in Ref. [H]. Here, we explain the TOF timing measurement scheme with the Time 
Stretcher [15] (TS) and TDC (TS-TDC), as it is essential to determining the bucket 
number of a colliding bunch for each event. 

The TOF readout system records a set of charges Qi and timings Tj from the rising 
edges of discriminator outputs for each PMT signal from the TOF detector. Figure lAll 
shows the block diagram of the TS-TDC scheme. The TS circuit finds the first rising 
edge T2 of the TS reference (reduced RF) clocks following the rising edge Ti of the 
TOF signal, and expands the time interval T21 = T2 — Ti by a factor of approximately 
/ = 20 to obtain T3 (T3 — T2=T2i x f+a). These times are then read out with the 
Belle standard TDC with a 0.5 ns least significant bit (LSB), providing a 25 ps LSB 
resolution. We reconstruct the TOF timing Tts-twc from timings, Ti, T2 and T3, using 



KERB Beam Collision Stability observed with the Belle Detector 



14 



the parameters / and a calibrated for each TS channeL We further apply a Time Walk 
Correction (TWC) to correct for timing variation proportional to 1/y/Q [T^]. 



Collision time 



Ref Tbco 



RF signal 
508.9 MHz 

Reduced RF signal 

BC/O reference time 
128 MHz 

TOF signal tjop 



Time Stretcher output 

to TDC1877S : LSB = 0.5 ns 



Common stop 

2.2\is - 30\is : given by GDL 




T3 T4 



Tstop 



Figure Al. Time Stretcher TDC scheme for TOF. The TS reference clock of 
approximately 8 ns period is generated from the KEKB RF signal of 508.9 MHz. 



The TS reference clock of approximately 8 ns (127.2 MHz) is generated from the 
RF clock of 508.9 MHz. As the RF clock is used for the entire KEKB timing control 
of the LER and HER beams, the TS clock is necessarily synchronized with the beam 
collision timing. 

Each event is triggered with a DAQ clock of approximately 16 ns period 
synchronized with every 8-th RF bucket. An RF bucket number wdaq (a multiple 
integer of 8) corresponding to the trigger is then recorded. In the standard TOF offline 
analysis, we determine a precise time relative to the DAQ clock. In particular, we 
determine an integer nxoF (a fine bucket number relative to nDAq); by taking into 
account track momenta and their path lengths to the TOF counter (times-of-fiight with 
e,/i, TT, kaon and proton mass hypotheses). We obtain nxoF with a candidate set of 
the mass assignments for each track, by optimizing a sum of squared time differences 

{STijf for all tracks in the event, where 6Tij = (TTs-TWCij - '^tof x Trf) - ^cxpij 
is the time difference using the j'-th mass hypothesis for the i-th track. Here, Trf is 
the RF period of approximately 2 ns. Thus, we find nxoF, and obtain the measured 
TOF timing tobsi = (^rs-TWCi ~ ''^tof x ^rf) for the track and the bucket number 
A'^buckct = '^DAQ + '^TOF of the collisiou bunch in the train. Figure [3] shows the number 
of yU-pair events as a function of the reconstructed bucket number. 

Appendix A. 2. Time resolution of the TOF detector 

Here, we discuss the sources of the time resolution for the TOF detector. Figure lA2] (a) 
shows the distribution of time difference (tA ~ ^b) between two measurements with 
independent channels tA and for the same input signal, in order to study the time 
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jitter in the TOF readout system. We estimate a time jitter of 35 ps from (T(tA — ^b)/ V^, 
which includes the time jitter in the discriminators of the TOF signal, errors in the TS 
(20 ps), and errors in the FASTBUS-TDC system. Figure lA2] (b) shows the distribution 
of time difference At = (tobs — ^exp) (=^ip) for /i-pair events. Here, t^hs is the measured 
time with TWC calibration after TS reconstruction, while texp is the predicted time 
calculated from the track trajectory in the CDC and its extrapolation to the TOF. 

Whenever both KEKB operation and the Belle detector (TOF and CDC) are stable, 
the mean of At should be zero. However, At varies run-by-run and in the long term. 
In practice, we correct for the run-by-run variation to set the mean At to zero for each 
run, so that we have no bias in physics analysis. 

The present time resolution of the TOF system is found to be 110 ps for yU-pair 
tracks. This resolution includes a systematic error of 40 ps consisting of the following 
components; a time jitter of 10 ps between the RF and TS clocks, 35 ps in the TOF 
readout system (front-end and TS-TDC circuits), 10 ps uncertainty in the time-of-flight 
calculation from /^-track information, and 13 ps due to the collision position spread of 
4.0 mm because of bunch length. The 10 ps uncertainty in the time-of-flight calculation 
is estimated for /i-pair events (by Monte Carlo), while it is larger than 10 ps for hadron 
tracks due to effects from hadronic interactions. The resolution of 110 ps is degraded 
from the initial resolution of 96 ps obtained in 1999 [H]. Taking into account the 
uncertainties discussed above, we extract intrinsic time resolutions of the Belle TOF 
detector of 102 ps and 87 ps for June of 2008 and 1999, respectively. The degradation is 
mainly due to ageing, a reduction of the attenuation length and light yield in the TOF 
scintillation counters over a ten year period. 




Figure A2. (a) Distribution of {t^ — ^b) for times measured by two independent 
readout channels A and B for the same input signal (b) Distribution of A^tof = 
^obs — ^oxp for /i-pair events for Exp 65, which gives a TOF time resolution of a= 
110 ps. 
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Appendix B. Beam bunch phase vs bucket number at KEKB 

Appendix B.l. Bunch-by-bunch phase detection 

In the KEKB collider about 450 beam position monitors (BPM) are mounted along the 
beam pipe for both of the LER and HER, in order to provide the beam position at 
specific locations. Each BPM consists of four electro-static electrodes to detect beam 
passages. One BPM is used to study possible beam phase variation as a function of 
bucket number in a beam train. The beam pick-up signal is fed into a gated beam- 
position monitor system (GBPM) that employs an I/Q (In-phase and Quadrature-phase) 
demodulator operating at the RF frequency of 508.9MHz to extract two orthogonal 
signals, Vsm and Kos for each bunch. A gate module in the front end selects a specific 
bucket in a beam train and enables bunch-by-bunch detection of the beam. The two 
orthogonal components and a specific bucket number are sampled by a beam revolution 
clock of approximately 10 /^s and stored in memory by GBPM. The beam phase (pheam 
of a specific bunch, with respect to its bucket phase 0rf synchronous to the 508.9 MHz 
RF clock, is calculated from the two components Vsin and Kos by: 

0boam - 0RF = tan^^(-^). (B.l) 

^^cos 

When the GBPM data are averaged over 2,000 turns, we estimate a resolution of ±0.1 
degrees in the phase measurement including the jitter of the RF reference clock. This 
corresponds to a time resolution of about ±0.5 ps. 



Appendix B.2. Transient beam loading effect 

The amplitude and phase of the RF in the accelerating cavities is modulated by the 
BAG, since the beam loading is different between the BAG and the beam train occupied 
by bunches. As a result, the synchronous position is shifted bunch-to-bunch in a beam 
train. Assuming that the cavity is operated at the optimum tuning and that the filling 
time of the cavity Tf is much longer than the revolution period Tq, i.e. T^/Tf ^ 1, the 
beam-phase difference A0ht between the head and tail bunches is given approximately 
in ref. [16] as: 

A</.HT = ^^UKFAtg, (B.2) 

where It is the total beam current, Rg is the shunt impedance, Vc is the cavity voltage, 
Q is the Q- value of the cavity, wrf is the angular RF frequency and Atg is a transit time 
of one BAG. In KEKB, the LER operates using ARES cavities with Tq/Tj = 0.54, while 
the HER operates using two different types of cavities, ARES and SC (about a half of 
the ARES Tq/Tj value). This causes the LER to have stronger non-linear behavior of 
the beam phase. 

In order to evaluate the modulation due to the BAG, a time domain simulation 
code was developed, taking into account beam loading on the accelerating mode of 
the cavities [17]. The simulation result was compared with the LER beam phase 
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measurement, as shown in Fig. IBll of Ref. [TB] . The simulation (dashed hne) accurately 
reproduces the data (closed circles) except for the first bucket. The rapid change in the 
beam head in the data is thought to be due to low Q longitudinal wakes that are not 
taken into account in the simulation. 
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Figure Bl. Beam phase (in degrees) as a function of bucket number. The sohd 
circles are the measured values while the dashed line shows the results of simulation. 
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